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Abstract 

This paper extends our recently developed Life Space Foam (LSF) model of moti- 
vated cognitive dynamics [1]. LSF uses adaptive path integrals to generate Lewinian 
force-fields on smooth manifolds, in order to characterize the dynamics of individ- 
ual goal-directed action. According to explanatory theories growing in acceptance in 
cognitive neuroscience, one of the key properties of this dynamics, capable of linking 
it to microscopic- level cortical neurodynamics, is its meta-stability and the result- 
ing phase transitions. Our extended LSF model incorporates the notion of phase 
transitions and complements it with embedded geometrical chaos. To describe this 
LSF phase transition, a general path-integral is used, along the corresponding LSF 
topology change. As a result, our extended LSF model is able to rigorously rep- 
resent co-action by two or more actors in the common LSF~manifold. The model 
yields substantial qualitative differences in geometrical properties between bilateral 
and multi-lateral co-action due to intrinsic chaotic coupling between n actors when 
n > 3. 

Keywords: cognitive dynamics, adaptive path integrals, phase transitions, chaos, 
topology change, human joint action, function approximation 



1 Introduction 



General stochastic dynamics, developed in a framework of Feynman path inte- 
grals, have recently [1] been applied to Lewinian field-theoretic psychodynam- 
ics [2,3,4], resulting in the development of a new concept of life-space foam 
(LSF) as a natural medium for motivational (MD) and cognitive (CD) psycho- 
dynamics. According to the LSF-formalism, the classic Lewinian life space can 
be macroscopically represented as a smooth manifold with steady force-fields 
and behavioral paths, while at the microscopic level it is more realistically rep- 
resented as a collection of wildly fluctuating force-fields, (loco)motion paths 
and local geometries (and topologies with holes). 
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A set of least-action principles is used to model the smoothness of global, 
macro-level LSF paths, fields and geometry, according to the following pre- 
scription. The action psycho-physical dimensions of Energy x Time ~ 
Effort and depending on macroscopic paths, fields and geometries (com- 
monly denoted by an abstract field symbol is defined as a temporal integral 
from the initial time instant tini to the final time instant tfin, 

sm = m dt, (1) 

with Lagrangian density given by 

£[$] = 

where the integral is taken over all n coordinates = x^{t) of the LSF, 
and are time and space partial derivatives of the variables over 

coordinates. The standard least action principle 

= 0, (2) 

gives, in the form of the so-called Euler-Lagrangian equations, a shortest 
(loco)motion path, an extreme force-field, and a life-space geometry of mini- 
mal curvature (and without holes). In this way, we effectively derive a unique 
globally smooth transition map 

F : INTENTIONt,^, ACTION,^^^, (3) 

performed at a macroscopic (global) time-level from some initial time tini to 
the final time tfin- In this way, we have obtained macro-objects in the global 
LSF: a single path described by Newtonian-like equation of motion, a single 
force-field described by Maxwellian-like field equations, and a single obstacle- 
free Riemannian geometry (with global topology without holes) . 

To model the corresponding local, micro-level LSF structures of rapidly fluc- 
tuating MD & CD, an adaptive path integral is formulated, defining a multi- 
phase and multi-path (multi-field and multi-geometry) transition amplitude 
from the state of Intention to the state of Action, 

{Action\I ntention) total V[w^e'^^^\ (4) 

^-^ LSF 

where the Lebesgue integration is performed over all continuous = paths+ 
fields + geometries, while summation is performed over all discrete processes 
and regional topologies The symbolic differential ©[w^] in the general 
path integral (22), represents an adaptive path measure, defined as a weighted 
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product 

N 

T>[w^] = lim Yl ^sd^l, (i = 1, n = con + dis). (5) 

°° s=l 

The adaptive path integral (4)-(5) represents an oo— dimensional neural net- 
work, with weights w updating by the general rule [1] 

new value{t -\- 1) — oldvalue{t) + innovation{t) . 



On the other hand, it is well-known that phase transitions (PTs) are phenom- 
ena which bring about qualitative physical changes at the macroscopic level 
in presence of the same microscopic forces acting among the constituents of a 
system. Their mathematical description requires to translate into quantitative 
terms the mentioned qualitative changes. The standard way of doing this is 
to consider how the values of thermodynamic observables, obtained in labora- 
tory experiments, vary with temperature, or volume, or an external field, and 
then to associate the experimentally observed discontinuities at a PT to the 
appearance of some kind of singularity entailing a loss of analyticity [22]. De- 
spite the smoothness of the statistical measures, after the Yang-Lee theorem 
[23] we know that in the N ^ oo limit non-analytic behaviors of thermody- 
namic functions are possible whenever the analyticity radius in the complex 
fugacity plane shrinks to zero, because this entails the loss of uniform con- 
vergence in N (number of degrees of freedom) of any sequence of real- valued 
thermodynamic functions, and all this depends on the distribution of the ze- 
ros of the grand canonical partition function. Also the other developments of 
the rigorous theory of PTs (see, e.g., [24,25]), identify PTs with the loss of 
analyticity. 

Similarly, experimental findings and theoretical insights of modern neuro- 
science converge on interpreting brain physiology within the conceptual frame- 
work of nonlinear dynamics, operating at the brink of criticality, which is 
achieved and maintained by self-organization [5]. In this approach, dynamical 
patterning of both brain activity and corresponding behaviors is examined in 
order to develop models of how brain and behavioral events are coordinated. 
Growing evidence supports the assumption that the linkages between events at 
microscopic level of neuronal assemblies and those at macroscopic behavioral 
levels are better explained as based on shared dynamics - not on any ontolog- 
ical priority [6,7]. This dynamics is characterized by met a-st ability and phase 
transitions in self-organized criticality [8] , both at the level of neuronal assem- 
blies [9,5], functional connectivity of the human brain [10] and corresponding 
behavior patterns [6,11]. A key feature of this approach is that phenomenolog- 
ical laws at the behavioral level can be connected to a field-theoretical descrip- 
tion of cortical dynamics [6] . Dynamic Field Theory (DFT) [12] extends this 
approach by developing field-theoretic representations of both behavior and 
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its environment [13], thus building on a long established, albeit metaphorical, 
tradition of behavioral force-field analysis [2,3,4]. Our LSF-formalism can be 
seen as a further extension of DFT. 

Regarding brain modelling: classical physics has provided a strong foundation 
for understanding brain function through measuring brain activity, modelling 
the functional connectivity of networks of neurons with algebraic matrices, 
and modelling the dynamics of neurons and neural populations with sets of 
coupled differential equations [14,15]. Various tools from classical physics en- 
abled recognition and documentation of aspects of the physical states of the 
brain; the structures and dynamics of neurons, the operations of membranes 
and organelles that generate and channel electric currents; and the molecu- 
lar and ionic carriers that implement the neural machineries of electrogenesis 
and learning. They support description of brain functions at several levels of 
complexity through measuring neural activity in the brains of animal and hu- 
man subjects engaged in behavioral exchanges with their environments. One 
of the key properties of brain dynamics are the coordinated oscillations of 
populations of neurons that change rapidly in concert with changes in the 
environment [16]. 

Also, most experimental neurobiologists and neural theorists have focused on 
sensorimotor functions and their adaptations through various forms of learn- 
ing and memory. Reliance has been placed on measurements of the rates and 
intervals of trains of action potentials of small numbers of neurons that are 
tuned to perceptual invariances and modelling neural interactions with dis- 
crete networks of simulated neurons. These and related studies have given a 
vivid picture of the cortex as a mosaic of modules, each of which performs a 
sensory or motor function; they have not given a picture of comparable clarity 
of the integration of modules (see [16] and references therein). 

The EEG analysis performed on rabbits and cats trained to discriminate con- 
ditioned stimuli in the various modalities (with EEG-recordings collected from 
high-density electrode arrays fixed on the epidural surfaces of primary sensory 
and limbic areas) has shown that cortical activity does not change contin- 
uously with time but by multiple spatial patterns in sequences during each 
perceptual action that resemble cinematographic frames on multiple screens 
[16]. The carrier waves of the patterned activity in frames have come in at 
least two ranges identified with beta (12-30 Hz) and gamma (30-80 Hz) oscil- 
lations. The abrupt change in dynamical state with each new frame, proposed 
to be formed by a phase transition [17], has not been describable either with 
classic integro-differential equations, or the algebras of neural networks. The 
initiation and maintenance of shared oscillations by this phase transition re- 
quires rapid communication among neurons. Several alternative mechanisms 
have been proposed as the agency for widespread synchrony. These are based 
in the dendritic loop current as the chief agent for intracellular communication 
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and the axonal action potential as the chief agent for intercellular communi- 
cation [16]. 

According to [16], many-body quantum field theory appears to be the only 
existing theoretical tool capable to explain the dynamic origin of long-range 
correlations, their rapid and efficient formation and dissolution, their interim 
stability in ground states, the multiplicity of coexisting and possibly non- 
interfering ground states, their degree of ordering, and their rich textures 
relating to sensory and motor facets of behaviors. It is historical fact that 
many-body quantum field theory has been devised and constructed in past 
decades exactly to understand features like ordered pattern formation and 
phase transitions in condensed matter physics that could not be understood 
in classical physics, similar to those in the brain. 

Communication by propagating action potentials imposes distance-dependent 
delays in the onset of re-synchronization during a phase transition over an 
area of cortex. The delays are measurable as brief but distance-dependent 
phase lags at the various frequencies of oscillation [17]. However, the length 
of most axons in cortex is a small fraction of observed distances of long-range 
correlation, with the requirement for synaptic renewal at each successive relay. 
These long-range correlations are maintained despite continuous variations in 
transmission frequencies that are apparent in aperiodic 'chaotic' oscillations. 

Some researchers have sought to explain zero-lag correlations with processes 
other than axo-dendritic synaptic transmission, stating that both electric fields 
and magnetic fields accompany neural loop currents. However, the electric 
potential gradients of the EEG have been shown by [18] to be inadequate in 
vivo to account for the long-range of the observed coherent activity, largely 
owing to the shunting action of glia that reduce the fraction of extracellular 
dendritic current penetrating adjacent neurons and minimize ephaptic cross- 
talk among cortical neurons. 

In this paper, to describe the LSF-phase-transitions, with embedded chaos, we 
use our adaptive path-integral (22) along the corresponding LSF-topology- 
change: 

(phase out ] phase in) := 4^ ©[tw^] e^'^'^l 

1/ topology — change 

This paper extends the earher establishment of the LSF model by introducing 
the study of chaos and phase transitions within this framework. This develop- 
ment is motivated - as was the original LSF model - by the potential for an 
improved theoretical basis for studying brain physiology, and thereby also for 
overcoming shortfalls in current artificial neural network function representa- 
tion and approximation technologies. 
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2 Geometrical Chaos and Topological Phase Transitions 



In this section we extend the LSF-formahsm to incorporate geometrical chaos 
and associated topological phase transitions. 

It is well-known that on the basis of the ergodic hypothesis, statistical me- 
chanics describes the physics of many-degrees of freedom systems by replacing 
time averages of the relevant observables with ensemble averages. Therefore, 
instead of using statistical ensembles, we can investigate the Hamiltonian (mi- 
croscopic) dynamics of a system undergoing a phase transition. The reason for 
tackling dynamics is twofold. First, there are observables, like Lyapunov expo- 
nents, that are intrinsically dynamical. Second, the geometrization of Hamil- 
tonian dynamics in terms of Riemannian geometry provides new observables 
and, in general, an interesting framework to investigate the phenomenon of 
phase transitions [21,35]. The geometrical formulation of the dynamics of con- 
servative systems [26] was first used by [27] in his studies on the dynamical 
foundations of statistical mechanics and subsequently became a standard tool 
to study abstract systems in ergodic theory. 

The simplest, mcchanical-like LSF-action in the individual's LSF-manifold E 
has a Riemannian locomotion form [1] 

1 f^fin ■ ■ 

^[q] ~ n ["^ij ^^i'' ~ ^(?)] (summation convention is assumed) 

(6) 

where aij is the 'material' metric tensor that generates the total 'kinetic en- 
ergy' of cognitive (loco)motions defined by their configuration coordinates 
and velocities q\ with the motivational potential energy V{q) and the standard 
Hamiltonian 

N 1 

H{p,q) = j:-p^ + V{q), (7) 
where Pi are the canonical (loco) motion momenta. 

Dynamics of N DOF mechanical-like systems with action (6) and Hamiltonian 
(7) are commonly given by the set of geodesic equations [32,33] 

^ + r;,^^ = o. (8) 

ds^ ^'^ ds ds 

where F*^ are the Christoffel symbols of the affine Levi-Civita connection of 
the Riemannian LSF-manifold E. 

Alternatively, a description of the extrema of the Hamilton's action (6) can 
be obtained using the Eisenhart metric [28] on an enlarged LSF space-time 
manifold (given by {(f = t,q^, . . . , q^} plus one real coordinate q^~^^), whose 
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arc-length is 



ds'^ = -2V{q){dqy + ttijdq'dq^ + 2dq^dq^^\ (9) 

The manifold has a Lorentzian structure [35] and the dynamical trajectories 
are those geodesies satisfying the condition ds^ = Cdt'^, where C is a positive 
constant. In this geometrical framework, the instability of the trajectories 
is the instability of the geodesies, and it is completely determined by the 
curvature properties of the LSF-manifold E according to the Jacobi equation 
of geodesic deviation [32,33] 

D'J\r.^ dq^ ,dq"^ _ 

^ + ^^,.^^^^-0, (10) 

whose solution J, usually called Jacobi variation field, locally measures the 
distance between nearby geodesies; D/ds stands for the covariant derivative 
along a geodesic and R j^m ^-^^ the components of the Riemann curvature 
tensor of the LSF-manifold E. 

Using the Eisenhart metric (9), the relevant part of the Jacobi equation (10) 
is given by the tangent dynamics equation [29,21] 

^ + itlVo^'^ = 0, (i = l,...,7V), (11) 

where the only non- vanishing components of the curvature tensor of the LSF- 
manifold E are 

R\ko = d^V/dq'dqK 

The tangent dynamics equation (11) is commonly used to define Lyapunov ex- 
ponents in dynamical systems given by the Riemannian action (6) and Hamil- 
tonian (7), using the formula [30] 

Ai = hm l/2Uog(Ef=,[J^(t) + JKt)\/^Ujm + m\)- (12) 

t— >oo 

Lyapunov exponents measure the strength of dynamical chaos. 

Now, to relate these results to topological phase transitions within the LSF- 
manifold E, recall that any two high-dimensional manifolds E^, and E^,/ have 
the same topology if they can be continuously and diffcrcntiably deformed into 
one another, that is if they are diffeomorphic. Thus by topology change the 
'loss of diffeomorphicity is meant [35]. In this respect, the so-called topological 
theorem [22] says that non-analyticity is the 'shadow' of a more fundamental 
phenomenon occurring in the system's configuration manifold (in our case the 
LSF-manifold): a topology change within the family of equipotential hyper- 
surfaces 

E, = {(g\...,g^)eR^in?\...,?^)=4, 
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where V and are the microscopic interaction potential and coordinates re- 
spectively. This topological approach to PTs stems from the numerical study 
of the dynamical counterpart of phase transitions, and precisely from the ob- 
servation of discontinuous or cuspy patterns displayed by the largest Lya- 
punov exponent Ai at the transition energy [30]. Lyapunov exponents can- 
not be measured in laboratory experiments, at variance with thermodynamic 
observables, thus, being genuine dynamical observables they are only be es- 
timated in numerical simulations of the microscopic dynamics. If there are 
critical points of V in configuration space, that is points Qc = [g^, . . . , g^] such 
that Vl^(g)|g^g = 0, according to the Morse Lemma [31], in the neighbor- 
hood of any critical point qc there always exists a coordinate system q{t) = 
[q^(t), .., q^{t)] for which 

Viq) = V{q,) -ql ~ql + fg+i + " " " + 9^, (13) 

where k is the index of the critical point, i.e., the number of negative eigenval- 
ues of the Hessian of the potential energy V . In the neighborhood of a critical 
point of the LSF-manifold S, (13) yields 

d^V/dq'dq^ = ±6,^, 

which gives k unstable directions which contribute to the exponential growth 
of the norm of the tangent vector J [30]. 

This means that the strength of dynamical chaos within the individual's LSF- 
manifold S, measured by the largest Lyapunov exponent Ai given by (12), is 
affected by the existence of critical points qc of the potential energy V{q). 
However, as V{q) is bounded below, it is a good Morse function, with no 
vanishing eigenvalues of its Hessian matrix. According to Morse theory [31], 
the existence of critical points of V is associated with topology changes of the 
hypersurfaces {S^j^gK- 

More precisely, let V^iqi, ■ ■ ■ ,qN) '■ R, he a smooth, bounded from 

below, finite-range and confining potential^]- Denote by S^, = V^^{v), v E R, 
its level sets, or equipotential hypersurfaces, in the LSF-manifold S. Then let 
V = v/N he the potential energy per degree of freedom. If there exists A^o, and 
if for any pair of values v and v' belonging to a given interval = [vq, Vi] and 
for any N > Nq then the sequence of the Helmoltz free energies {F/v(/3)}ArgN 
- where P = 1/T (T is the temperature) and (3 E Ip = {(5[vq) , (3{yi)) - 
is uniformly convergent at least in C'^{Ip) [the space of twice different iable 
functions in the interval //j], so that limjv_^oo -^Af ^ C^i^fs) neither first 
nor second order phase transitions can occur in the (inverse) temperature 



These requirements for V are fulfilled by standard interatomic and intermolecular 
interaction potentials, as well as by classical spin potentials. 



8 



interval {P{vo), P{vi)), where the inverse temperature is defined as [22,35] 

(3{v) =dS'^j^\v)/dv, while s'^j^\v) = N-Hog [ d^q 

Jv{q)<vN 

is one of the possible definitions of the microcanonical configurational entropy. 
The intensive variable v has been introduced to ease the comparison between 
quantities computed at different A?"- values. 

This theorem means that a topology change of the {E^j^gK at some Vc is a 
necessary condition for a phase transition to take place at the corresponding 
energy value. The topology changes implied here are those described within 
the framework of Morse theory through 'attachment of handles' [31] to the 
LSF-manifold S. 

In the LSF path-integral language [1], we can say that suitable topology 
changes of equipotcntial submanifolds of the individual's LSF-manifold S can 
entail t her mo dynamic-like phase transitions [37,38,39], according to the gen- 
eral formula: 

(phase out | phase in) := ? P[w$] e''^'*'. 

•^topology— change 

The statistical behavior of the LSF-(loco)motion system (6) with the standard 
Hamiltonian (7) is encompassed, in the canonical ensemble, by its partition 
function, given by the phase-space path integral 



Zn = £p_,,^b]^[?] exp |i \pq - Hip, q)] drX , (14) 



where we have used the shorthand notation 

dq{T)dp{T) 



Jtoi 



V\p]V[q] ^ I n 



top-ch J 27r 

The phase-space path integral (14) can be calculated as the partition function 
[36], 



ZnW) = / n dp^dq^e-f'^'M = M c^g^e-^^(«) 

1=1 Vr^/ i=l 

/ \ — 



where the last term is written using the so-called co-area formula [20] , and v 
labels the equipotcntial hypersurfaces E„ of the LSF-manifold E, 

E„ = {(g\...,g^)eR^|l-(g\...,g^)=t;}. 
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Equation (15) shows that the relevant statistical information is contained in 
the canonical configurational partition function 

Note that Z'^ is decomposed, in the last term of (15), into an infinite summa- 
tion of geometric integrals, 

/ da/WVVl 

defined on the Once the microscopic interaction potential V{q) is 

given, the configuration space of the system is automatically foliated into 
the family {Si,}„g]R of these equipotential hypersurfaces. Now, from standard 
statistical mechanical arguments we know that, at any given value of the 
inverse temperature (3, the larger the number A^, the closer to = 
are the microstates that significantly contribute to the averages, computed 
through Zn{P), of thermodynamic observables. The hypersurface is the 
one associated with 

Uf, = {Z^)-' ll[dq^V{q)e-f'''^'^\ 

the average potential energy computed at a given /3. Thus, at any /3, if is 
very large the effective support of the canonical measure shrinks very close 
to a single S„ = S„^. Hence, the basic origin of a phase transition lies in a 
suitable topology change of the {S^}, occurring at some Vc [36]. This topology 
change induces the singular behavior of the thermodynamic observables at a 
phase transition. It is conjectured that the counterpart of a phase transition 
is a breaking of diffeomorphicity among the surfaces S^, it is appropriate to 
choose a diffeomorphism invariant to probe if and how the topology of the S^, 
chang function of v. Fortunately, such a topological invariant exists, the 

Euler characteristic of the LSF-manifold S, defined by [32,33] 

N 

X(S) = E(-1)'&^(S), (16) 

k=0 

where the Betti numbers bkC^) are diffeomorphism invariants!^ This homolog- 
ical formula can be simplified by the use of the Gauss-Bonnet-Hopf theorem, 
that relates x(^) with the total Gauss-Kronecker curvature Kg of the LSF- 
manifold S 

X(S)= / Kg da, (17) 

where 

da = \J det{a)dx^dx^ ■ ■ ■ dx^ 

^ The Betti numbers are the dimensions of the de Rham's cohomology vector 
spaces H''{T,]R) (therefore the bk are integers). 
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is the invariant volume measure of the LSF-manifold S and a is the deter- 
minant of the LSF metric tensor aij. For technical details of this topological 
approach, see [34]. 

The domain of validity of the 'quantum' is not restricted to the microscopic 
world [19]. There are macroscopic features of classically behaving systems, 
which cannot be explained without recourse to the quantum dynamics. This 
field theoretic model leads to the view of the phase transition as a conden- 
sation that is comparable to the formation of fog and rain drops from water 
vapor, and that might serve to model both the gamma and beta phase transi- 
tions. According to such a model, the production of activity with long-range 
correlation in the brain takes place through the mechanism of spontaneous 
breakdown of symmetry (SBS), which has for decades been shown to describe 
long-range correlation in condensed matter physics. The adoption of such a 
field theoretic approach enables modelling of the whole cerebral hemisphere 
and its hierarchy of components down to the atomic level as a fully integrated 
macroscopic quantum system, namely as a macroscopic system which is a 
quantum system not in the trivial sense that it is made, like all existing mat- 
ter, by quantum components such as atoms and molecules, but in the sense 
that some of its macroscopic properties can best be described with recourse 
to quantum dynamics (see [16] and references therein). 

Phase transitions can also be associated with autonomous robot competence 
levels, as informal specifications of desired classes of behaviors for robots over 
all environments they will encounter, as described by Brooks' subsumption 
architecture approach [45,46,47]. The distributed network of augmented finite- 
state machines can exist in different phases or modalities of their state-space 
variables, which determine the systems intrinsic behavior. The phase transition 
represented by this approach is triggered by either internal (a set-point) or 
external (a command) control stimuli, such as a command to transition from 
a sleep mode to awake mode, or walking to running. 



3 Modelling Human Joint Action 

Cognitive neuroscience investigations, including flVIRI studies of human co- 
action, suggest that cognitive and neural processes supporting co-action in- 
clude joint attention, action observation, task sharing, and action coordination 
[40,41,42,43]. For example, when two actors are given a joint control task (e.g., 
tracking a moving target on screen) and potentially conflicting controls (e.g., 
one person in charge of acceleration, the other - deceleration), their joint per- 
formance depends on how well they can anticipate each other's actions. In 
particular, better coordination is achieved when individuals receive real-time 
feedback about the timing of each other's actions [43] . 
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To model the dynamics of the two-actor joint action, we propose to asso- 
ciate each of the actors with an n— dimensional Riemannian LSF-manifold 
S, that is a set of their own time dependent trajectories, = {a;*(tj)} and 
= {y-'(tj)}, respectively. Their associated tangent bundles contain their 
individual nD (loco)motion velocities, TS^ = = dx^/dti} and TS^ = 

{y-^{tj) = dy^/dtj}. Further, following the general formalism of [1], outlined 
in the introduction, we use the modelling machinery consisting of: (i) Adap- 
tive joint action at the top-master level, describing the externally-appearing 
deterministic, continuous and smooth dynamics, and (ii) Corresponding adap- 
tive path integral (22) at the bottom-slave level, describing a wildly fluctuat- 
ing dynamics including both continuous trajectories and Markov chains. This 
lower-level joint dynamics can be further discretized into a partition function 
of the corresponding statistical dynamics. 

In particular, by extending and adapting classical Wheeler-Feynman action- 
at-a-distance electrodynamics [44] and applying it to human co-action, we 
propose a two-term joint action: 



^[^^ y'l ^i; ij\ 



with 



x\U)-y^{t, 



a^(3j6{Itj) x\ti)y\tj) dUdtj + - j^gijx\t)x\t) dt 

2 



where IN <ti,tj,t < OUT. 



(18) 



The first term in (18) represents potential energy of the cognitive /motivational 
interaction between the two agents and (3j\2l It is a double integral over a 
delta function of the square of interval P between two points on the paths in 
their Life-Spaces; thus, interaction occurs only when this interval, representing 
the motivational cognitive distance between the two agents, vanishes. Note 
that the cognitive (loco)motions of the two agents aj[x*(tj)] and Pj[y^(tj)], 
generally occur at different times ti and tj unless tj = tj, when cognitive 
synchronization occurs. 

The second term in (18) represents kinetic energy of the physical interaction. 
Namely, when the cognitive synchronization in the first term takes place, the 
second term of physical kinetic energy is activated in the common manifold, 
which is one of the agents' Life Spaces, say = {x\ti)}. 

Conversely, if we have a need to represent coaction of three actors, say Oj, j3j 
and ■y^ (e.g., a, in charge of acceleration, jSj - deceleration and 7^— steering), 
we can associate each of them with an nD Riemannian Life-Space manifold, 
Sq, = {x\ti)}, = {y-^itj)}, and = {z'^(tfc)}, respectively, with the corre- 
sponding tangent bundles containing their individual (loco)motion velocities, 

^ Although, formally, this term contains cognitive velocities, it still represents 'po- 
tential energy' from the physical point of view. 
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T^a = {x'iU) = dxydU},TEf, = {y^{t,) = dy^/dt,} and TS^ = {z\tk) 
dz^ /dtk}. Then, instead of (18) we have 



A[U, tj, tk;t] = - / / ai{U)(3^ (t,) 7, (4)5(4^) x'{U) y^{tj) z\tk) dUdtjdtk 

+ 1 [ W.^it, 1, q) (f(l" dt, where IN < U, tj, h, t < OUT (19) 
2, Jt 

with 4, = [x\t;) - yKh)f + Wih) - + [z\t,) - x\u)]\ 

Due to an intrinsic chaotic couphng, the three-actor (or, n— actor, n > 3) 
joint action (19) has a considerably more comphcated geometrical structure 
then the bilateral co-action (18)o It actually happens in the common 3nD 
Finsler manifold Ej = Eq, U U E^, parameterized by the local joint coordi- 
nates dependent on the common time t. That is, Ej = {q^{t)^ r = 1, ...,3n}. 
Geometry of the joint manifold Ej is defined by the Finsler metric function 
ds = F{q^, dq^), defined by 

F'^il^ q) = drsil, (l)Q^(f 1 (where Qrs is the Riemann metric tensor) (20) 

and the Finsler tensor Crst{q,q), defined by (see [32,33]) 

_ l d^F^{q,q) _ 1 dgrs , . 

^rst[q,q) ^g-rgqsQqt 2dq-dq^- ^ ^ 

From the Finsler definitions (20)-(21), it follows that the partial interaction 
manifolds, Eq, U E^, E^ U Ej^ and Eq, U E^, have Riemannian structures with 
the corresponding interaction kinetic energies, 

1 . . 1 . 1 . 

Tap = ^9iji'y\ Tay = -gikx'z^, Tp^ = -gjky^z^. 



At the slave level, the adaptive path integral (see [1]), representing an 00— dimen- 
sional neural network, corresponding to the adaptive bilateral joint action (18), 
reads 

{OUT\IN) := fv[w, X, y] e^^t^'S^^*-*^], (22) 

where the Lebesgue integration is performed over all continuous paths = 
x*(tj) and y^ = y-'{tj), while summation is performed over all associated dis- 
crete Markov fluctuations and jumps. The symbolic differential in the path 
integral (22) represents an adaptive path measure, defined as a weighted prod- 
uct 

N 

V[w, X, y] = \im Y[wl:dx'dy\ (i, j = 1, n). (23) 

s=l 

^ Recall that the necessary condition for chaos in continuous temporal or spatio- 
temporal systems is to have three variables with nonlinear couplings between them. 
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Similarly, in case of the triple joint action, the adaptive path integral reads, 

{OUT\IN):= fv[w;x,y,z;q]e'^^''^'^^'''''\ (24) 
with the adaptive path measure defined by 

N 

V[w;x,y,z;q] = lim Y[ Wi^f.^dx'dyHz^dq\ {ij,k = l,...,n; r = l,...,3n). 

(25) 

The adaptive path integrals (22) and (24) incorporate the local Bernstein 
adaptation process [48,49] according to Bernstein's discriminator concept 

desired state SW {t + 1) = current state IW{t) + adjustment step AW (t). 



4 Discussion 

This paper has developed an adaptive path integral approach to modelling 
topological phase transition, chaos and joint action in the LSF-manifold. The 
traditional neural networks approaches are known for their classes of func- 
tions they can represent This limitation has been attributed to their low- 
dimensionality (the largest neural networks are limited to the order of 10^ 
dimensions [53]). The proposed path integral approach represents a new fam- 
ily of function-representation methods, which potentially offers a basis for a 
fundamentally more expansive solution. 

This new family of function- representation methods is now capable of repre- 
senting input/output behavior of more than one actor. However, as we add the 
second and subsequent actors to the model, the requirements for the rigorous 
geometrical representations of their respective LSFs become nontrivial. For a 
single actor or a two-actor co-action the Riemannian geometry was sufficient, 
but it becomes insufficient for modelling the n-actor (with n > 3) joint action, 
due to an intrinsic chaotic coupling between the individual actors' LSFs. To 
model an n-actor joint LSF, we have to use the Finsler geometry, which is a 
generalization of the Riemannian one. This progression may seem trivial, both 
from standard psychological point of view, and from computational point of 
view, but it is not trivial from the geometrical perspective. 

The robustness of biological motor control systems in handling excess de- 
grees of freedom has been attributed to a combination of tight hierarchical 
central planning and multiple levels of sensory feedback self-regulation that 

^ Here we are talking about functions in an extensional rather than merely inten- 
sional sense; that is, function can be read as input/output behavior [54,55,56,57]. 
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are relatively autonomous in their operation [50]. These two processes are con- 
nected through a top-down process of action script delegation and bottom-up 
emergency escalation mechanisms. There is a complex interplay between the 
continuous sensory feedback and motion/action planning to achieve effective 
operation in uncertain environments, such as movement on uneven terrain 
cluttered with obstacles. 

Complementing Bernstein's motor control principles is Brooks' concept of 
computational subsumption architectures [45,47], which provides a method 
for structuring reactive systems from the bottom up using layered sets of be- 
haviors. Each layer implements a particular goal of the agent, which subsumes 
that of the underlying layers. Similar architectures have been proposed to ac- 
count for the mechanism of cognitive and motor working memory, sequence 
learning and performance (see, e.g., [11]). According to [45,47], a robot's low- 
est layer could be "avoid an object" , on top of it would be the layer "wander 
around" , which in turn lies under "explore the world" . The top layer in such 
could represent the ultimate goal of "creating a map". In this config- 
uration, the lowest layers can work as reflexive mechanisms, while the higher 
layers contain control logic implementing more abstract goals. 

The substrate for this architecture comprises a network of finite state ma- 
chines augmented with timing elements. A subsumption compiler compiles 
augmented finite state machine descriptions into a special-purpose scheduler 
to simulate parallelism and a set of finite state machine simulation routines. 
The resulting networked behavior function can be described as: 

final state w{t + l) = current state w{t) + adjustment behavior f{ Aw {t)). 

The Bernstein weights, or Brooks nodes, = in (23) are updated 

by the Bernstein loop during the joint transition process, according to one of 
the two standard neural learning schemes, in which the micro-time level is 
traversed in discrete steps: 

(1) A self-organized, unsupervised (e.g., Hebbian-like [51]) learning rule: 

<{t + 1) = Kit) + ^(wtfit) - wtf{t)\ (26) 

where a = (j{t), rj = ri{t) denote signal and noise, respectively, while new 
superscripts d and a denote desired and achieved micro-states, respec- 
tively; or 

(2) A certain form of a supervised gradient descent learning: 

wtj{t + l) = wt^{t)-r)VJ{t), (27) 
where 77 is a small constant, called the step size, or the learning rate, and 
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VJ(n) denotes the gradient of the 'performance hyper-surface' at the 
t—th iteration, 

where t = to, ti, ...,ts then t + 1 = ti,t2, t^+i. 

Both Hebbian and supervised learninj^ are used in local decision making 
processes, e.g., at the intention formation phase (see [1]). Overall, the model 
presents a set of formalisms to represent time-critical aspects of collective 
performance in tactical teams. Its applications include hypotheses genera- 
tion for real and virtual experiments on team performance, both in human 
teams (e.g., emergency crews) and hybrid human-machine teams (e.g., human- 
robotic crews). It is of particular value to the latter, as the increasing auton- 
omy of robotic platforms poses non-trivial challenges, not only for the design 
of their operator interfaces, but also for the design of the teams themselves 
and their concept of operations. 

Some specific problems that Brooks poses in [47] include: (i) how to combine 
many (e.g., more than a dozen) behavior generating modules in a way which 
lets them be productive and cooperative; (ii) how to automate the building 
of interaction interfaces between behavior generating modules, so that larger 
(and hence more competent) systems can be built; and (iii) how to automate 
the construction of individual behavior generating modules, or even to au- 
tomate their modification. An aggregation of phase-transition-related order 
parameters within individual actors need to be triggered (either by an internal 
or external control mechanism) into collective alignment to be able to per- 
form the joint action. The 'guiding' forces help to guide the better alignment 
or fine-tunning of individual LSFs for useful co-action outcomes to emerge. 
Using the combined LSF-geometry approach proposed in this paper should 
provide some useful answers to the above questions. 

Here we remark that collective phase transitions, which have more coupled de- 
grees of freedom than individual ones, necessitate more stringent constraints, 
to avoid/control a higher-dimensional chaos. The sophisticated chaos-control 
techniques, including constraining contextual boundaries - choosing a target 
subspace of the joint LSF-manifold - and guiding force-fields, need to be de- 
fined to alow desired collective behaviors to emerge from both chaotic and 
non-chaotic sets of possible initial evolution alternatives. 



^ Note that we could also use a reward-based, reinforcement learning rule [52], in 
which system learns its optimal policy: 

innovation{t) = \reward{t) — penalty {t)\. 
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5 Conclusion 



Extending the LSF model to incorporate the notions of meta-stabihty, phase 
transitions and embedded geometrical chaos has enabled representation of in- 
creased complexity in goal-directed action, including the capability to repre- 
sent joint action by two or more co-actors. This capability remains consistent 
with modern neuroscience theorizing linking macro-behavioral meta-stability 
and phase transitions to microscopic-level cortical neurodynamics. 

There is a degree of correspondence between phase transition mechanisms 
for cognitive performance in humans and transitions between stable behavior 
states and competency levels in autonomous robots. The approach developed 
in this paper offers a theoretical framework to integrate observations and mod- 
els of both individual and collective robot behaviors and competencies, capable 
of coping with the increased complexity of the real world. This framework can 
both guide future substantive empirical work into collective robot behaviors 
and be vahdated by it. 

The new model developed in this paper offers substantial improvements over 
the geometrical properties regarding multiple-actor systems, due to the chaotic 
coupling between the actors. We have also discussed how the proposed path 
integral represents a new family of function representation techniques that 
may expand on the range of function types that are currently afforded by 
standard neural-network models. Specifically, we are interested here in bio- 
logically plausible function representation as a means for characterizing the 
input/output behavior of multi-actor systems, and thus regard function rep- 
resentation in an extensional sense. The full realisation of these possibilities 
in practical applications is a subject of ongoing research. 
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